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ABSTRACT 


A theory  is  developed  for  the  parametric  excitation  of  coupled  waves 
in  a uniform  magneto  plasma.  Instead  of  using  electrostatic  approximation 
which  is  valid  only  for  longitudinal  waves,  a general  formulation  of  the 
parametric  coupling  equations  is  developed  by  using  Hamiltonian  approach. 
The  plasma  is  pumped  by  a monochromatic  electromagnetic  wave  with 
frequency  which  decays  into  two  decay  waves  and  (e.  g.  , a 
plasmon  and  a phonon)  with  Either  one  of  the  plasma  modes 

coupled  with  the  pump  wave  (through  the  linearized  Vlasov  equation)  to 
induce  a polarization  which  may  act  as  the  source  of  another  mode.  A 
general  procedure  to  calculate  the  induced  polarizations  is  introduced  by 
transforming  the  linearized  Vlasov  equation  into  the  oscillating  frame  of 
reference,  in  which  the  equation  can  then  be  solved  by  the  method  of 
characteristics.  The  threshold  power,  initial  growth  rate  and  the  frequency 
shift  (due  to  finite  initial  growth  rate)  of  the  parametric  instabilities  can  be 
obtained  from  the  coupled  mode  equations  and  the  results  obtained  can  be 
applied  to  all  the  wavelength  region.  Several  interesting  cases  are  pre- 
sented and  compared  favorably  with  experiments. 

The  simultaneous  excitation  of  second  harmonic  of  electron- cyclotron 
and  electrostatic  ion-cyclotron  oscillations  by  means  of  a monochromatic 
electromagnetic  wave  at  <*>^“9.  23  GHz  has  been  observed  experimentally. 
The  experiment  has  been  conducted  in  a microwave  sustained  beam  plasma. 
By  switching  the  pump  power  between  two  levels,  the  time  evolution  of  the 
excited  electrostatic  ion-cyclotron  wave  is  studied.  In  order  to  explain  the 
observed  nonlinear  saturated  state  of  the  instability,  the  nonlinear  damping 
rate  due  to  two  possible  mechanisms  has  been  examined.  Harmonic  gener- 
ation gives  very  small  contributions,  therefore,  anomalous  diffusion  caused 


by  the  coherent  oscillation  of  the  electrostatic  ion-cyclotron  wave  plays  an 
important  role.  Working  in  the  guiding  center  frame,  the  effective  non- 
linear damping  frequency  pertaining  to  the  diffusion  process  is  obtained. 
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I.  INTRODUCTION 


Plasma  parametric  instabilities  occur  when  the  pump  amplitude 

exceeds  a threshold  which  depends  on  the  set  of  coupled  modes  and  on 

their  damping  effects.  It  has  become  of  special  interest  following  the 

recognition  that  such  processes  play  important  roles  in  the  absorption 

of  electromagnetic  power^.  Experimental  results  showed  that  parametric 

instabilities  may  form  a channel  for  efficient  transfer  of  electromagnetic 

energy  into  a rather  large  size  hot  plasma,  and  the  nonlinear  effects  of 

these  instabilities  may  introduce  a mechanism  which  heats  plasma  more 

2 

efficiently  thc.n  classical  collisional  absorption  . Therefore,  most  of  the 
experimental  effort  has  been  expended  on  observing  the  plasma  heating 
rates.  Besides  plasma  heating,  theoretical  and  experimental  studies  of 
parametric  decay  instabilities  will  also  yield  information  about  the  linear 
and  nonlinear  damping  rates  of  the  instabilities,  and  the  development  of 

3 

plasma  turbulence.  In  the  previous  report  , a detailed  experimental  study 
of  the  dynamics  of  the  process  had  been  given.  Klein  and  Cheo  ’ introduce 
a technique  for  experimentally  determining  growth  and  decay  characteristics 
of  parametrically  excited  plasma  waves.  They  found  that  the  dominate 
linear  damping  mechanism  of  the  ion  acoustic  wave  in  their  plasma  is  due 
to  ionization  collisions,  not  due  to  ion- neutral  collisions  which  is  the  usual 
phenomenon.  Theirs  is  a unmagnetized  plasma,  while  this  effort  is  the 
continuation  of  their  experimental  study  into  a magnetoplasma.  Therefore, 
different  modes  are  involved.  Since  there  are  no  available  theories  for  our 
experiment,  a general  theory  is  also  developed  by  using  system  total 
Hamiltonian  approach  to  deal  with  the  problem.  Comparison  with  our  ex- 
periment has  also  been  made. 

Indeed,  the  detailed  ways  of  implementation  of  the  ideas  of  parametric 

2 

decay  instabilities  to  plasma  heating  is  still  developing.  Grek  and  Porkolab  , 
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Porkolab  et  al.  , and  Okabayashi  et  al.  had  performed  experiments  by 

shining  microwave  power  onto  a plasma  column  in  the  extraordinary  mode 

of  propagation  in  the  region  of  the  upper  hybrid  frequency.  Wavelength 

measurements  of  the  decay  waves  showed  decay  into  upper  hybrid  waves 

(and  Bernstein  waves)  and  lower  hybrid  and/or  ion  acoustic  waves.  It  was 

shown  that  significant  ion  and  electron  heating  occurred  only  above  threshold 

6>  T 8 

for  parametric  instabilities.  There  were  also  several  experiments”’  ’ 
performed  in  the  regime  of  lower  hybrid  frequency  which  showed  parametric 
decay  into  lower  hybrid  waves  and  ion  acoustic  waves,  ion  cyclotron  waves, 
ion  quasi  modes,  or  drift  waves  may  occur.  Again,  strong  plasma  heating 
was  observed  above  threshold  for  parametric  instabilities  (both  ion  and 
electrons).  Experimental  technique  to  observe  the  plasma  heating  rates  is 
by  using  a multigrid  energy  analyzer  to  measure  the  modified  electron  and 
ion  distribution  function.  Another  technique  is  also  developing  by  switching 
the  pump  power  between  two  levels  and  using  the  Fabry-Perot  interferometer 
to  measure  the  ion  temperature  growth  rate. 

In  the  parametric  process,  the  incident  pump  wave  is  couples  to  two  or 
more  natural  modes  of  the  plasma  by  satisfying  both  conservation  of  frequency 
and  conservation  of  wave  vector  relations.  The  frequency  mixing  allows 
energy  supplied  to  the  system  at  one  frequency  to  be  converted  to  another. 

The  physical  origin  of  mode  coupling  mechanism  is  fairly  easy  to  understand 
on  a qualitative  basis.  From  linear  wave  analyses,  the  dielectric  tensor  of 
the  plasma  is  density  dependent.  Hence  it  will  be  modulated  by  the  presence 
of  the  vibrational  modes  of  the  plasma,  resulting  in  a mixing  action  such  that 
the  input  radiation  and  vibrational  modes  beat  together  to  produce  sum  and 
difference  frequencies.  These  in  turn  may  act  as  sources  for  parametric 
excitations  if  the  frequency  and  wave  vector  matching  conditions 
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are  satisfied. 

By  assuming  exp{-iu)t)  time  behavior  for  the  mode  of  interest,  decay 
or  growth  of  the  mode  will  be  decided  by  the  dispersion  relation  of  the 
system.  Therefore,  most  of  the  theoretical  investigations  are  directed 
toward  finding  the  modified  dispersion  relations  with  the  presence  of  a 

q 

pump  wave  in  the  system.  The  first  wo^k  is  done  by  Silin.  His  theory  is 
largely  based  on  the  hydrodynamic  equations  for  a cold  plasma.  DuBois  and 
Goldman  analyze  the  parametric  coupling  of  Langmuir  and  ion-acoustic 
oscillations  based  on  a Green's  function  perturbative  method  (harmonic 
approximation),  which  is  restricted  to  the  case  when  the  radiation-induced 
energy  of  the  particles  is  small  compared  to  their  thermal  energy.  They 
showed  that  the  plasma  can  be  unstable  to  certain  applied  frequencies 
for  pump  power  above  a threshold,  a characteristic  of  parametric  excitation. 
Later,  Nishikawa,  and  Lee  and  Su^^  used  hydrodynamic  model  and  obtained 
same  conclusions.  However,  in  order  to  justify  the  hsrmonic  approximation 
for  greater  intensities,  Jackson  pointed  out  that  it  is  important  to  estimate 
the  range  of  frequencies  which  produces  instabilities  when  the  intensity  is 
large. 

Parametric  coupling  of  electrostatic  waves  in  a magnetized  plasma 

have  been  treated  by  Aliev  et  al.  by  Amano  and  Okamoto,^^  and  by 
16  IT 

Porkolab,  ’ and  others.  These  authors  used  the  linearized  Vlasov 
equation  with  self  consistent  potential  field,  which  can  be  solved  easily 
after  transformation  to  the  oscillating  frame  of  reference.  Alieve  et  al. 
found  the  frequency  range  of  a pump  wave  in  which  the  plasma  is  unstable 
turns  out  to  be  much  broader  than  in  the  case  of  a unmagnetized  plasma. 
Amano  and  Okamoto  extend  the  theory  to  inhomogeneous  cases.  Porkolab 
analyzed  the  resulting  dispersion  relation  of  Aliev  et  al.  in  the  limit  of 
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weak  coupling,  and  obtained  threshold  powers  to  excite  the  upper  and 
lower  hybrid  (or  ion  acoustic)  modes  simultaneously.  He  also  introduced 
a new  type  of  kinetic -dissipative  instability.  All  of  the  above  investigations 
are  based  on  electrostatic  approximation  for  the  plasma  modes  and  hence 
valid  only  for  longitudinal  waves.  In  the  presence  of  a static  magnetic 
field,  most  of  the  plasma  modes  become  hybrid,  and  the  electrostatic  approxi- 
mation holds  only  in  the  limiting  case. 

In  this  investigation  a general  formulation  of  the  parametric  coupling 

equations  in  a homogeneous  magnetized  plasma  is  developed  by  using 

1 8 

Hamiltonian  approach.  The  coupling  coefficients  of  the  parametric 
equations  are  derived  from  the  collisionless  Boltzman-Vlasov  equation. 

With  transformation  method  and  trajectory  techniques,  the  induced 
polarization  currents  are  derived.  Applications  to  various  forms  of 
modal  coupling  are  given,  and  comparison  with  experimental  works  is 
made  whenever  possible. 

A set  of  experiments  in  the  Poly  HCD  plasma  was  performed.  A 
microwave  horn  is  mounted  at  the  end  of  the  plasma  beam.  In  order  to 
obtain  a quiescent  plasma,  a new  technique  to  create  an  ultrastable  micro- 
wave  sustained  plasma  is  introduced.  When  the  pump  signal,  at  9.  23  GHz, 
is  above  the  threshold  level,  the  second  harmonics  of  the  electron  cyclotron 
wave  and  the  electrostatic  ion-cyclotron  wave  are  simultaneously  excited. 

The  time  evolution  process  of  the  electrostatic  ion-cyclotron  wave  is 
studied  by  switching  the  microwave  pump  power  from  one  level  to  another. 
Linear  damping  rate,  initial  growth  rate,  and  the  threshold  power  can  be 
obtained. 

In  addition, the  nonlinear  saturation  mechanisms  observed  in  the 
experiment  are  discussed.  An  analysis  by  using  the  guiding  center 
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technique  was  made  and  it  is  established  the  nonlinear  damping  is  due 
to  the  anomalous  diffusion  produced  by  the  ponderamotive  forces  of  the 

L’. 

instabilities.  Comparison  with  experiments  is  favorable. 
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n.  GENERAL  FORMULATION  FOR  THE  PARAMETRIC  EXCITATION 


1.  Overall  Physical  Features 

It  is  well  known  that  the  vibrational  modes  of  a medium  will  modulate 

the  permittivity  or  polarizability  of  the  medium.  Assume  that  there  exist 

two  normal  modes  or  waves  ,S  and  L (e.  g.  , a phonon  and  a plasmon)  in 

the  medium  with  frequencies  and  in  the  linear  regime;  then  they 

will  modulate  the  polarizability  of  the  medium  resulting  from  the  density 

perturbation  of  the  medium.  The  problem  can  be  treated  mathematically 

by  expanding  the  electrostrictive  polarizability  in  a Taylor's  series  in 

terms  of  normal  mode  vibrational  field  5E  related  to  S and  6E.  related 

— s —i 

to  L 


.a-ij  =afij+(g^  5Es^ 

SO" 


where  i,  j,  o-  =1,2,3  or  x,  y,  z 


where  the  Einstein  summation  convention  is  understood.  Hence  the  total 
induced  polarization  due  to  the  presence  of  a third  wave  field  E is 


= 1 E =P%6P..6/». 

Z (T 


5^1=  (a^^o  ®^s<r^j  ^^e  the  induced 

scr  •'  i a J 

polarization  associated  with  the  Raman  and  Brillouin  effects.  Here  we 

have  assumed  IeI  >>  ISE^  | and  |^}  >>  |SE^  | 


2.  Hamiltonian  Approach 

We  first  normalize  both  S and  L waves.  Assume  the  particles 
associated  with  S and  L have  reduced  masses  M and  m , momentum 
K and  JC,  displacements  Q and  U respectively.  The  displacement  and 


momentum  variables  are  the  canonical  conjugate  pairs.  Hence  the  xinper- 

(k)  (k) 

turbed  Hamiltonian  density  H and  H','  are 

•'os  oi 


H<^^  4 S [-2— 

OS  Z (j  ’■  M 


+ Mco'  (k)Q^{k)] 


H 


(k)  1 


oi 


Z?  t 


,2 


m 


+ mu^ 


(2.  1) 


(2.2) 


I I 

where  Wg(k)  (}S)  resonant  frequencies  of  the  system -without 

the  pump  field,  and  the  field  associated  with  the  corresponding  modes  are 


6E  (k)  = - — ^ a)^(k)  Q(k) 

K 


n, 

— 

e Jn~ 
o 


5E,{k)  =f 


(k)  U(k)  + — ^ Mk)  x^z 


(2.  3) 


(2.  4) 


where  Wg(k.)  (k)  the  resonant  frequencies  with  the  presence  of 

pump  field,  n^  is  the  background  electron  (or  ion)  density.  and 
are  cyclotron  frequencies  of  the  two  species  including  the  sign  of  the 
charge.  Here  we  also  assume  that  a uniform  magnetic  field  is 

applied,  and  the  plasma  is  singly  ionized.  It  is  easy  to  show  that  the 
Poisson  bracket  relations  hold: 


(a,&  . 

K^(k))=  5(k,k)  6^^ 

(2.  5) 

= 5<k,k')  6,^ 

(2.  6) 

The  pump  field  is  assumed  to  be  spatially  homogeneous  (dipole  approxima- 
tion). Hence  the  induced  polarizations  corresponding  to  the  frequency 
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components  in  question  can  be  expressed  as 


— ) 6E^  (k)  eT  +. 

dSE^  (k)  ° 3 

io-  — 


(k)  8s,..  . ^ 

5P  =( il— ) 6E^^(k)  E,  + c.c. 

^ 86E"  (k)  ° J 

SC  “ 


(2.7) 


(2.  8) 


where  the  superscripts  denote  respectively  exp  (+iwt).  Then  the  inter- 
action Hamiltonian  densities  for  S and  L are 


^'s^^  =^Pi^^5E^.(k)  = 6p(^^^  [M  wJ(k)Q.(k)  + (K(k)xz)J 


e Jn 
^ o 


(2.  9) 


= - SP^^^^SE^  .(k)=-  2^yU.(k)+^^  (if(k)xz)J 


e Jn~ 

^ o 


(2.  10) 


The  total  Hamiltonian  densities  of  the  S and  L are 


pj(k)  _ ^(k)  ^ fj'(k) 

S OSS 


njk)  ^ j^(k)  ^ (k) 

i oi  i 


(2.  11) 


(2.  12) 


It  is  known  in  general  that  the  equation  of  motion  for  any  cononical  variable 
A is 


[d_  + 2 r ^ + r ^ 1 A = { {A,  H},  H} 
dt 


(2.  13) 


where  1 / r is  the  phenomenological  relaxation  time,  and  H is  the  total 
Hamiltonian  density.  (2.  13)  together  with  (2.  1)  to  (2.  12)  can  be  used  to 
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derive  the  coupled  mode  equations; 

Evaluation  of  the  Poisson  bracket  { {Q^(k),  with  the 

aid  of  (2.  1),  (2.  5),  (2.  9)  and  (2.  11)  yields 


{{Q.{k), 


hW}. 
s ^ 


'2  (k) 

.,^(k)Qi(k)--^6pW 

e Jn 

^ o 


Therefore  the  appropriate  equation  of  motion  for  Qj(k)  is  given  by  (2.  13), 

+2r  Q.(k)  + {r  (k))  Q.(k)  = - 5P'^*?'  (2.  14) 

1 3 1 9 S 1 I"  SI 

e JrT 

with  the  dot  denoting  . 

Similarly  by  evaluating  { {K.(k),  with  the  aid  of  (2.  1), 

(2.  5),  (2.  9)  and  (2.  11)  we  have; 


{ {K  (k)  , } = - «/(k)K.(k)  - Ma>;,^(k)  — ^ (1  X 5pW) 

1 9 9 S*  I •■*5  1 


Q. 


(k), 


e,/n 


and 


K.(k)+2r^K.(l^+(r^+w'^^(k)  )K.(k)  =-  — (z  x (2.  15) 

e Jn~ 


Q. 


Equation  of  motion  in  the  variable  SE^(k)  can  be  obtained  by  combining 
(2.  14)  and  (2.  15)  with  the  aid  of  the  relation  (2.  3).  The  result  is 


dt  n e 


Mw^(k)  Mw',^(k)nf 

® 5P'*''+ 


— s 


-■--■y  - (lx  5p[^^^x£ 


n e 
o 


(2.  16) 


Similarly  for  the  L mode  we  obtain; 
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and  combining  the  relevant  components,  equations  (2.  18)  and  (2.  19)  can  be 
rewritten  as 


^2  A 2 Mw'^(k)  Mw'^(k)n? 

[■S+2r  ^ + « (k)]6E:(k)  = ^_6p-(k)  + (g  X 6^(k)  )x2 
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d“  d 2 + nriw.  (k) 


<2  2 

_ , mw  (k)n  , 

= — S-  5P.(k)+ ^-^(«x5P,(k))xg 

dt  n e n e 

° ° (2.21) 


Here  (2.  20)  and  (2.  21)  fit  the  standard  form  of  parametric  coupled  wave 
eqviations  with  6P  remaining  formal,  which  is  to  be  derived  in  the  next 
section. 


m.  KINETIC  APPROACH  TO  CALCULATE  THE 
COUPLING  COEFFICIENTS 


1.  Transformation  of  the  Vlasov  Equations  to  the  Oscillating  Frame  of 
Reference 


To  derive  the  coupling  coefficients  ) and  for 


'a5E,  o 
io- 


'asE  'o 

sc 


determing  5P,  we  start  from  Vlasov  equation  for  a homogeneous  back- 


ground distribution  function  f 


(o). 


a 


e , 


a 


a 


8v 


(3.  1) 


or 


where  e^,  m^  and  v^  are  respectively  the  charge,  mass,  and  velocity  of 


a particle  of  species  q,.  The  corresponding  solution  can  be  written  in  the 
fo  rm 


f^°\v,t)=f  (v  - ^ 7^  dt' R (t-t')  • E(t'),t)  =f  (v,t) 
a a a°  a —a  — a°~ 


(3.2) 


where 


R (t)  = 

~a 


cos  n t 
a 


sin  Q t 
a 


sin  O t 
a 


cos  Q t 0 

a 


(3.3) 


® <3o 
« me 


v = V - — “ 
— —am 


r ' 

J dt*  R (t-t')  • E(t') 


(3.  4) 


Hence  in  the  oscillating  frame  equation  (3.  1)  reduces  to 


® 1 

aO  , a 1 


+r:f^  4vx  B,  • —02,: =0 


at 


me—  1 
a 


dv 


(3.  5) 


ijs/tsaumumsmixii 


".r:' 


J 


-13- 


Since  any  function  of  |v|  satisfies  (3.  5),  therefore,  let's  assume 
f (v,t)  = f (v,,,v,  ),  a Maxwellian  distribution. 

<yO  ~ aO  — 

It  can  be  shown  easily  that  to  the  first  order  equation  for  the  pertur- 
bative distribution  function  6f  (v  ,t)  is 

a ~'a  a 


aSf  (v  , r , t) 
Q-n  a 

at 


+v 

~a 


aSf  {v  , r ,t)  e , 

— 2: — + _2LrE(t)+lv  xB  ] 

ar  m '•  — ' ' c —Q  — o ' 

~a  at 


aSf  (v  , r ,t) 
a ~a 2 

9v 

"“a 


9f^°\v  ,t) 

a zo 

dv 

~a 


= 0 


(3.6) 


For  convenience  we  transform  the  original  variables  (v  ,£  , t)  to  (v,  r,  t) 

a 01 


of  the  oscillating  reference  frame  and  introduce  the  function 


ijj  (v,  r,t)  = 5f  (v  , r ,t)  = 6f  (v+-^  dt'R  (t-f)  • E(t'),  r + f , t) 
a a a a o-  “ ~ ~a 

(3.7) 

where  t_  =jr+i 

Of  a 

and  i dt"  S (t'-t")  • E(t")  (3.8) 

.»  -a  ~ 

with 


g aSf  (v  ,r  ,t)  e 

. -tt-g  -ai-  /*  dl'  B (t-f)) 


asf  (v  , r ,t) 

a ~o  ~a 

dv 

~ot 


+ 


asf  (v  ,£  ,t) 

— a Tg g 

ar 

~a 


From  (3.  6),  the  equation  for  d*  (v,  r,  t)  is  then 

a 
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_a_ 

ar 


+ n (vx4)  • -^1  4>^(v,  r,  t)  = •> 


a 


a 


^ 5E(r+i  ,t) 


af  ^(v) 

av 

(3.  9) 


2.  Integration  along  the  unperturbed  trajectory  in  the  phase  space. 

Equation  (3.  9)  can  now  be  solved  by  the  method  of  characteristics. 

In  this  method  the  perturbative  distribution  function  is  calculated  in  the 
Lagrangian  system  of  coordinates.  Let's  define  an  unperturbed  trajectory 
of  the  particles  as 


dr 

dF 


= V 


dv 

-5^  =n  (vxz) 
dt  a 


The  solution  of  (3.  10)  is 


(3.  10) 


v'  = R (t'-t)  • V (3.  11) 

— —a  ~ 

and 

r'  = r+^  L (t'-t)  • V (3.12) 

— — —at  — 

a 


where 

L (t)  = 

~a 


and 


r = r(t) 


sin  Q t 
a 

-(1-cosn  t) 
a 

0 


1 - cos  SI  t 
a 

sinSl  t 
a 

0 


0 

0 


n t 
a 


= n r^R  (T)dT 

or  o 


(3.  13) 


r'  = r(t')  , v=v(t)  , v;'  =v(t') 


Hence  along  the  unperturbed  trajectory  (3,  10)  equation  (3.  9)  may  be 
written 


^ <3. 14) 


-•3-^ 
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Substituting  (3.  7)  in  (3.  14),  we  obtain 


j e 

^ 5f  (v+-f-  dt-  a (t-f)  • E(t')  , r + i ,t)  = 
at  Of  — m ■'  —a  — — —a 

a "“c 


e 9f  (v) 

Of 


(3.  15) 


Because  of  the  spatial  homogeneity  of  the  background,  we  may  assume  a 

ik*  jr 

coordinate  dependent  of  the  non- equilibrium  increment  6f  — e and 

ik  •_£ 

the  perturbed  field  5^  - e Then  (3.  15)  becomes 


(v  + 


O' 


ik-  {r+l  ) 

dt'R  (t-t')- E(t'),k,t)e  “ 

—Of  ~ 


e 

a 


df  (v) 

Qi£ 

9v 


ik • (r+f  ) 


(3.  16) 


Let  us,  however,  restrict  ourselves  to  unstable  (growing)  solutions.  We 
can  thus  invert  this  differentiation  by  integrating  (3.  16)  with  respect  to 
time  along  the  unperturbed  particle  trajectory  as 


dt'E  (t-f)  .E(t'),k,t)  = - 
~ ~ 


5E(k,t') 

-oc 


df 


dv' 


ik-(r'.r)  %[i,(t')-i^(t)] 


df 


(3.  17) 


From  equation  (3.  12),  let 

(t)  = k • (r  - r ')  = k • L (-t)  • v and  t = t - t' 

or  — — — — — (y  — 

a 


' . ..aaaw.- 


(3.  18) 
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Then  equation  (3.  17)  becomes 


6f  (v  -I — ^ 
a “ m 


^ .T-'/j.ix  1_  i.\  — Ot  p ^ 


r dt's  (t-f)  •E(t'),k,t)  = 

~Q 


■SC 


m 


a o 


af 

• { fcos  T + 6)x  + sin  (n  t + 0)  y ] 

oV|  a a 


9f  -i4,  (t)  ik.[i  (t-T)-i  (t)] 

+ ^ S } e " e ® dr 

dv„ 


where  we  have  used  following  relations 


V = V,  cos  0 X + V,  sin  0 y + v,,  z 


and 


v'  = V,  cos  (Q  T + 0)  3c  + V,  sin  ( il  t + 0)  y + Vn 

^ /V  /V 


then 


9fjv')  8f  ^(v) 

ja^o 


^ [ cos  (n  T + e)  X + sin  (il  T + 0)  y ] +■ 


8f 


8v'  av,  ' a a 


8v 


S!° 


Let  the  time  dependence  of  the  pump  electric  field  be 


E(t)  = 2E  , sin  u t + 2E-,  cos  u t 
— ' ’ — 1 o — 2 o 


then 


k-ff  (t-T)  -1  (t)  ] =-V  f*  dt'  /*'  dt-'k-E  (t'-t")-  E(t") 


a 


t-T  -OC 


= (a  , +b  ,)  [ sin  0)  t - sin  co  ^(t-T)] 

01  ^ Qt^  **0  O 

+ (a  ,-b  ,)rcosu  t - cos  CO  ^(t-T)l 

' Ot^  ® ^ 


(3.  19) 


(3.  20) 


(3.21) 


/ 

i 
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^ dt'E  (t-t')-E{t')  =i:^  _J!o  [A  . (ET+e!)-  a*.  (E^  + E2)  1 
m *'  — Q rn  ^ ••  — ct  — a ^ ^ 

Qi  -oc  a w " ^ 

° “ (3.22) 


where 


e k • A • E.  e f? 

a . = 2 — ^ ^ , b . = 2 ^ 5-  '5;  X k • E.  , i = 1,  2 

/yi  m 2 ai  m ,2  ^ — — 1 

“ /yW-n  /y(i)((0-n) 

o a o o a 


, 0)  - i2 

_ =xx+“yy-i z z 

— ^rv  ^ 

Cl) 

O 


= i e 


^^2  "— 2 ® 


-ico  t 
o 


, io)  t 

eJ  = - i E j e ° 


, ioj  t 

^2=^26  ° 


2 ^2 

S2  u)  - 

A =xSi:  + i-^(xy-yx)  + — z z 

— Of  (*> 

O 0) 

o 


A is  the  complex  conjugate  of  A . 
— a —a 


We  may  then  calculate  the  total  induced  current  from  (3.  19)  as 


5j(k,t)  = In  e /dv[v+— a:-  dt'R  (t-f) 

ofaror  ni  ■<  —a 

a ”* 


E(t')16f  (v+-^  /^dt'E  (t-t')*E(t'),k,t) 


af  — m ■'  —Q 
a -» 


(3.23) 


(3.24) 


(3.  25) 


(3.26) 
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CO  = U)  + w . 

oat 


for  the  conservation  of  energy 
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iSg  = - = is  for  the  conservation  of  momentum  (dipole  approximation) 

where  and  are  the  characteristic  frequencies  of  two  plasma  modes  in 
a collisionless  plasma  with  the  presence  of  pump  field. 

Therefore,  the  induced  current  density  for  the  S (or  L)  mode  due  to 

the  coupling  between  L (or  S)  mode  and  the  pump  wave  can  be  obtained  with 
the  aid  of  (3.  27)  and  (3.  28)  as 


k- A . (eT+e!) 

**  —/V  —1  — ^ 


^ 9f  o 9f 

/vdv/*dT{-^  [cos  (n  -H-e)5c+sin(n  -p^e)^  ]+_o° 

o a a ^ ' 9v,,  ^ 


OJ 


-iw.T  i(i,  (t)  iw  t ) 

e (1-e  ° )j  . 64(k.t) 


(3.  29) 


and 


+ in  e oj  Qf 

%(iS.t)=Z— ^ “2^  -(Ej+E^)/ dv /"°dT{-^[cos(n  T+e)Sc 

m w ur  i ^ Q, 


a o a 

9f  iw  T -id  (x) 

+ sin(n  Tfe)y  ] + -x5^'4}e  ® e ® 

a 9v,,  ‘ 


k- A*-  (Et+Et) 
. rr  ~g  — 1— ^ 

^ 0) 


9f 


9f 


I [cos(n  .T+e)jc+sin(n  T+e)trl+-^  ^ } 

o a a ^ ' dv  „ ^ 

-icj  T 

“ (1-e  °)}-6E;(k.t)  (3.30) 


"i<^_(T)  -ICJT 

e e 


In  carrying  out  those  integrations,  we  specifically  choose  the  wave  vector 
k on  the  x-z  plane  and  write 


is  " k^  S:  + k,,  'i 


(3.  3 1) 
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Equation  (3.  18)  becomes 


<p  (t)  [ sin(f2  T+  e)  - sin  6 1 k,,v,,T 

ot  a a 


(3.  32) 


where 


i = k,  V,  /fi 

— — Cl 


(3.  33) 


then  we  have 


+ 0C  +00 


exp  [-i(;)  (t)  ] - Y.  E J ( 4?  )J  )exp{-i[n(Q  T+e)-n'e+k„v,,T  1} 
“ n=-ocn'=-oc  » 


(3.  34) 


Therefore,  we  obtain  following  relation 

9f  df  .....  . 

f dT{5-2-  [ cos(nT+0)x+sin(nT+0)y]  +•5^  z}  e 


1 Z . 

n=-oc  n 


n*(p,Z^;n)[l-W(ZJ]  - 2 z} 


(3.  35) 


for  a Maxwellian  distribution  function  f 

o 

where 


k„(T/m)2 


, (3=^T/mf2  and  A^(p)  = l^(p)  e"P  (3.36) 


a(p,Zn;n)=i^A^((3)x3t+  [^A^(p)-2pA;^(p)  iry  + z;A^(p)V2-inA;((3){Sc5r-ySt) 
+ ^ Z^A^{p)(xH45t)+  Z^A;(p)(y^-‘iy)  (3.  3 7) 


Substituting  (3.  35)  in  (3.  29)  and  (3.  30),  the  results  are 
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3 ^ 

mT(w-n)  I'-nJ-i  a 


a a o a 


W(Z“  ) , 

+ i3  A' (3  )V1  + ' A„(3  U 

Kii  ^ '^Q  J flf 


3 

1.  or  O' 

V 

1 

k,  u 

J-  o 

Zj 
^ n 

k,  (it),  +nI2  ) 
-‘-I  a 

3 

1 ’^a  Of 

[7 

1 

Ufa  7*^  , 

k,  u 

J-  o 

n 

k,  (w  -nfl  ) 

■*•  8 a 

- ^a'  ^sn' 

(e; +E-) 


(3.  38) 


3 

new 

sj#  (k.t)=i:  — ^ ^-2°--. 

^ “ m T (w  -q"^) 

a a ^ a 


«E;(k,  t)-|[  I [ 1 -W(Z“„)l[„A„(p^)S 


I W(z®  ) . ^ 

+ i3  A„{p  )y] ir^A  (p  )z  A 

^Of  K,,  n ^cx  J —or 


0 Q^r  1 2:®  - * 

Z T— 7 — rr-;n(3  , z"  ;n)[l-W(Z°'  )]  - r— [k- A ] 
k.w  l„  k,  (w^-nn  ) - Pff'  sn’  sn'-*  k.w  {'■—  ==a 

J-  o n * s ^ s J 


^2 

■ - W(Z»„)1 

J-o'-nJ-'i  cx 


kj_u^  ] — -a  — 1 — 


(3.  39) 


where 


Z°  = 

i n 


u,+nn  u -nfi  , 

^ , -^1  . 3=k2T/mn 

k.,(T  /m  k,.(T  /m  a ^ c^J  a 

a/  a a/  a 


whence  the  induced  polarizations  are  obtained  as 


6P"(k,t)=i-6j'(k,t)  = 6^  (k,t)  • (— ^)  • E' 

— s-  Wg  -s-  -i  - 95E+ 


6^  (k.t)  = - 6J^  (k,t)  = 6E:(k.t)  . (— ^ ) • 

-I  - WjJ  -i  - S gg^- 


where 


_ _ - + + + 

E =E.  +E„  and  E =E.  +E, 
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IV.  DETAILED  ANALYSIS  OF  THRESHOLD  POWER  AND  . 
INITIAL  GROWTH  RATE 

1.  Conditions  and  Properties  of  the  Instabilities  Occurring  on  the  Basis 
of  a General  Form  of  the  Coupled  Mode  Equations 

In  the  following  paragraph  we  will  show  that  the  coupled  wave  equa- 
tions (2.  20)  and  (2.  21 ) derived  in  Section  II  can  be  reduced  to  the  following 
scalar  form  similar  to  those  introduced  by  Nishikawa^^: 


dt 


4.  or  ±4.  2 


X = i\YE‘(t) 


(4.  1) 


dt 


1 .u  or  ± ^ 2 

2 dt  “i 


Y = -i^XE^(t) 


(4.  2) 


where  for  simplicity,  we  treat  only  the  case  with 
Xu  = real  > 0 

Taking  the  Fourier  transform  of  (4.  1)  and  (4.  2),  we  obtain 
2 


r 2 

W - OJ 


+ 2ir  w 
s s 


X(w)  = XE^Y(cj-oj^) 


- wj  + 2ir^(u-w^)j  Y(w-w^)  = jiE^X(w) 


(4.  3) 


(4.4) 


where 


+00 


X(t).  Y(t)=  / |^e''‘^X(«).  Y(w). 

-00 

Setting  the  determinant  of  the  coefficient  matrix  equal  to  zero,  the  disper- 
sion relations  which  determine  the  frequency  and  damping  (or  growing)  of 
the  waves  under  consideration  are  obtained,  namely. 
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r + 2ir  (w-w  )^  - + 2irfl  (w-Gj  ) - XuE^  = 0 (4.  5) 

[s  sJ_o  £ ioj'^o 


If  the  solution  of  (4.  5)  is  written  in  the  form 


Gj  = X + iy 


(4.  6) 


then  X and  (-y  ) are  the  frequency  and  damping  rates,  respectively,  of  the 
new  normal  mode  with  parametric  coupling  effect.  It  becomes  unstable  if 


y > 0 


(4.7) 


Considering  first  the  threshold  case  in  which  y = 0,  = E^  and  (4.  5)  can 

be  separated  into  two  equations 


(X-GJ 


+r_g(x-u.J(x2-.J)=  0 


(X^-GJ^)  (x-w^)^  - - 4rgr^  x{x-g;^)  - = 0 


(4.  8) 


(4.  9) 


Since  F and  F,,  are  arbitrary,  the  threshold  intensity  of  the  pump  field 
s Jc 


and  the  frequencies  of  X and  Y are  given  by 


(4.  10) 


X - GJ  = -GO. 

o £ 


(4.  11) 


(4.  12) 


Thus  under  the  condition  of  frequency  match,  there  is  no  frequency  shift 
of  X and  Y at  the  threshold. 


With  the  pump  field  > E^ , we  no  longer  can  set  y = 0.  Sub- 
stituting (4.  6)  in  (4.  5)  and  separating  the  real  and  imaginary  parts,  we 


obtain 
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(x^-y^-w2_2p^yj  j^(x-u^)2  . a;2-2r^y  - 4x(x-w^){y+r  )(y+r^ 


- = 0 


(4.  13) 


x(y+r  ) (x-w^)2  . . Y^Y  + (x-u)^)(y+r^)(x2-y2-wj-2rgy) 


(4.  14) 


If  the  frequency  shift  is  small  compared  with  w , then  (4.  1 4)  can  be  ap- 


proximated as 


-2  Wg(x-Wg)+r^y  I to^(y+rg) 


2 2 2 _ -rr" -.'-i'j 

X -y  - -2r^y= , 0 


and  (4.  13)  becomes 


(4.  15) 


4(y+rg)(y+r_g  )t^g  = Xi^ej 


After  some  algebraic  manipulation,  we  find 


_ -(r3+r^)  V(r,-r  )^-HX^E2yc.^c.^ 


(r,^rpH-/(r^-r3)^  + 4r3r^Ef/E'f 


y{y+2r  ) 

X =:  W + jr 2- 

s 2u 

s 


(4.  16) 


(4.  17) 


(4.  18) 


In  general  is  much  larger  than  F^.  therefore  we  may  express  (4.  17) 
and  (4.  18)  in  the  vicinity  of  the  threshold  as 


y ^s  ^ ^ 

c 


(4.  19) 
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X a 


ui 

s 


UI 

s 


u 


s 


1 


(4.  20) 


The  initial  exponential  growth  rate  and  the  frequency  shift  are  linearly 

proportional  to  the  punap  intensity.  These  results  are,  however,  valid  only 

when  the  frequency  shift  is  small  compared  to  and  the  pump  intensity  is 

near  the  threshold.  These  phenomena  had  actually  been  observed  in  Stenzel 

20 

and  Wong's  experiment  in  an  unmagnetized  plasma. 


2.  Parametric  Decay  into  Longitudinal  Modes  (electrostatic  approximation) 
Thus  we  can  write 


6E”(k,  t)  = -ik0g{k,  t)  and  6^^(k,  t)  = 


(4.  21) 


and  substitute  these  relations  in  (2.  20)  and  (2.  21  ).  We  have 


d^  d 2 

Mw'^(k) 

0 (k,  t)  = 2 2-’ 

a(k*  a) 

72  ^ ^ 

+ ; 

dt 

nek 

o 

a 

O’ 

(k,  t) 


(k)n.^ 

s ~ -L  V 

2,  2 - 
nek 
o 


zx 


0a 


• E’  \ 0^(k.  t) 


(4.  22) 


mw^(k) 

0«  (Is*  - 2 2 — 

nek 

a(k  • a) 

86  E" 

o 

— s 

E''’0  (k,  t) 


mw'  ^ (k)ng 

^ 2,  2 - 
nek 
o 


z X 


Oor 


06  E 


^■^^■(k,  t)  (4.23) 


(a)  When  the  pump  wave  is  the  ordinary  mode  and  the  decay  waves 
are  ion  acoustic  mode  and  electron  plasma  mode 


Let  _E(t)  = 2E„  cos  t z (i.  e.  , = 0 


k=k.  (l.e.,  k||=k,  k_^=0.  0^=0) 
and  the  substitution  of  (4.  24)  into  (4.  22)  and  (4.  23)  will  yield 


(4.  24) 


2 9 f 

f?  + f + -J(k)  .1 » z W(- 


k/ T^/m^ 
a a 


)\  E'6^^k,  t) 


(4.  25) 


" 4 3 

+ 2r.  4 +-;.^(k)  <pUh.t)=  ) 


W(-  ^ )>Et(fe~(k.  t) 


(4.  26) 


Since  the  decay  modes  we  are  considering  are  in  the  frequency  domains 
such  that 


V T./M  « w /k  « ■/  T /m  « w. /k 

1 s e i 


then  (4.  25)  and  (4.  26)  become 


j 2 1 n \ kew 

72  ^ I <Pl{h <^;{k. t)E-{t) 

dt  I mcj 

J o 


(4.  27) 


(4.  28) 


ittWIt! 


^+^2  0;(k.t)  = 


k 


(4.  29) 


where  we  have  used  relations 

2 2 

9 9 9 9 4jm  e „ „ 4ffn  e 

= k"^T  /M,  = 2 — ^ J-  and  kj  = — (4.  30) 

s e i pe  m o d T ' ’ 

e 

Equations  (4.  28)  and  (4.  29)  had  been  obtained  by  Nishikawa^^,  who  worked 
with  a hydrodynamic  model.  If  we  take  the  Fourier  transform  of  (4.  28) 

and  (4.  29),  the  results  would  be  the  same  as  DuBois  and  Goldman^ 

11  12  13 

Nishikawa  , Lee  and  Su  , and  Jackson's  conclusions. 

Therefore,  the  threshold  field  is  given  by 


8 / Trn  T 
o e 


where  E.,  = 2E„. 

th  2 


u,  w. 
s i 


(b)  When  the  pump  field  is  extraordinary  mode  (Ej.B  ) and 

In  I » u r ° 

e pe 


In  this  case  we  may  write 


^ = 2E  (xcos  u t - y sinu)  t)  (i.  e. , E„  = xE  , = - E y) 

O O O O X o 


k = k X + k z 
X II 


Substituting  (4.  32)  in  (4.  22)  and  (4,  23),  we  obtain 


d2  d 2 

^ + 2r  77  + w,  <^J]S.t)  = (k  6 P’,  + k 6 P’  ) 

dt^  s dt  s ^s  - ^ ^2^2  ' j.  si  II  s3 

J ° 

iMw*  2^? 

+ l-TT^  k 6P', 


nek 
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2,2  X si 


(4.  31) 


(4.  32) 


(4.  33) 
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where 
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(4.  35) 
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(4.  37) 
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(4.  38) 


(i)  Decay  waves  are  electrostatic  ion  cyclotron  wave  (ion 
acoustic  wave)  and  upper-hybrid  wave 

Thus,  we  may  assume-^^  « k /k^  « 1 such  that  \V'(Z^^)=  0, 

and  W(Zg^)  = 0 for  all  n and  W(Z®^)  = 0 for  n + 0,  W-(Z®^)x\^{^^)  = 1. 

With  these  approximations,  equations  (4.  35)-(4.  3 8)  may  be  re- 


duced to 


0  

2 2 
2w  . mT  k j 

1 e d 


n e k 


5 (1  + k^/  k^) 

2mT  J ^ ^ 

e x 


(4.  39) 


and  (4,  33)  and  (4.  34)  become 
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2 JO  iMew  , 9 9 + 

^ A . ,, 0.,  t, 

e £ (4,40) 
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(4.  41) 
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Therefore,  the  threshold  field  can  be  calculated  as 


8V  7rn  T 
o e 


J „ OJ 

si  pe 


„ „ k'^T  /M 

where  E^.  = (2E  K.  , ^ 

th  ' oth'  s 1 

X a 

2 1 

had  been  obtained  by  Porkolah 


k^T  /M 

9 9-  and  k /kf 

1+kVk:  ^ 

I d 


(4.  42) 


a 1.  This 


(ii)  Decay  waves  are  lower  hybrid  and  upper  hybrid  waves 

In  this  case  k /k  «_/^.  therefore  W(Z  )=  0 for  all  n. 
11  X yM  n' 

\P.  \ fp  2 ~ J' 

Since  = -j-  +Y(-^  ) + ~ \nj  '^J^=  “i  +‘^s' 

k2  ^2  ^2 

Kj  16  6 
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CO  s:  ^ cj 

pe  ‘ e ' s 


and  = |(!j  + o^Aj(|3J  ; kf  „ ' 


With  these  approxinnations  , equations  (4.  35)-  (4.  38)  become 
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(4.  43) 


and  (4.  33)  and  (4.  34)  become 
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(4.  45) 


Therefore;  the  threshold  field  can  be  calculated  as 


th 
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^s“i  -/2w  fn 
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(4.  46) 


which  also  agrees  with  Porkolab's  result 


I 21 
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where  E.,  = (2E  ) , , w“  = il 
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(c)  For  a forced  linearly  polarized  pump  field  with  upper 
hybrid  and  lower  hybrid  waves  as  decay  waves 


In  this  case  we  write  E = 2E  cos  co  t x (i.  e.  E„  = xE  ; E,  = 0) 

— o o —2  o — 1 


and  k = kx+kz  (k/k  «-,/'r7  ) 

- i M II  1 M 


Since  ^ = k'^T  /mi2^  * « 1 and  w « w .. 


o e pe  e j. 


e pe  e 
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With  these  approximations;  equations  (4.  33)  and  (4.  34)  become 
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Therefore,  the  threshold  field  is  given  as 


8 V irn  T 
o e 


^s^Ji  VTu), 


(4.  49) 


3.  Parametric  Decay  into  Electrostatic  Ion  Cyclotron  Wave  and  Har- 
monics of  Electron  Cyclotron  Wave  (Hybrid  Mode)  in  a Uniform 
Magneto  Plasma 


For  the  harmonics  of  electron  cyclotron  modes,  |k|  is  small 

(|k|  «k,)  such  that  /3  « 1,  and  W(Z®  )«0  for  all  n. 
a e in 

Since  the  electrostatic  ion  cyclotron  mode  is  longitudinal  wave',  we 

may  assume  6E'(k,t)=  -ik0'(k,t),  where  k = k x + k z with 

Am  I Therefore,  W(Z®  ) ==  0 for  n ^ 0 and  W(Z®  ) A (/3  ) «,  1. 

V II  X sn  so  o e 

Because  the  nonlinear  coupling  due  to  electrons  of  the  plasma  is 
much  stronger  than  due  to  ions;  therefore  we  may  neglect  the  ion  terms 
of  equations  (3.  4 1)  and  (3.  42).  The  resulting  components  of  the  induced 
polarizations  can  be  simplified  to 
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(4,  50) 
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and  the  coupled  mode  equations  are  given  by  equations  (4.  33)  and  (2.  21).  From 
(4.  52),  (4.  53)  and  (2.  2 1),  we  obtain  the  relation  6ET„(k,  t)  = i—  dE"^,  (k.  t). 

X ^ ^ 00  ^ 1 '—  ' ' 

o 

thus  the  dispersion  relation  for  the  L mode  is  given  as 


, liVVe'.^e  , 


(4.  55) 


let  u)^  = Infill  (1  + <2i^)  where  n>  2 and  substitute  into  (4.55),  we  have 
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Since  the  polarization  of  the  harmonics  of  electron  cyclotron  wave  can  b« 


determined  from  linear  dielectric  tensor  as 
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thus  the  linearly  polarized  pump  field  can  be  written  as 
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Therefore,  the  resulting  coupled  mode  equations  are  given  as 
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and  the  threshold  field  can  be  calculated  as 
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4.  Comparison  Between  Theory  and  Experiments 


a.  In  the  case  of  paranaetric  decay  into  electrostatic  ion  cyclotron 
wave  and  second  harmonics  of  electron  cyclotron  wave  in  a uniform  mag- 
neto plasma.  (My  experiment) 

With  the  aid  of  (4.  60),  the  value  of  parameters  can  be  found  from 


Table  2,  the  threshold  field  is  calculated  as 


E , = 1 3.  7 volts/ cm 
tn 


Near  the  microwave  horn,  the  pump  field  is  almost  uniformly  distributed 
on  the  cross-section  of  the  chamber,  hence  the  theoretically  predicted 


threshold  pump  power  may  be  calculated  as 

^2 


Threshold  pump  power  a x Trr*^  x c =:  44  watts 

OTT 


where  r = 7.  5 cm  is  the  radius  of  the  chamber. 

The  threshold  pump  power,  experimentally  determined,  is  about 
47.  1 watts  as  shown  in  Figure  10.  This  value  is  very  close  to  the  theo- 
retically predicted  value,  and  the  extra  3 watts  is  believed  to  be  used  to 
sustain  the  plasma. 

b.  In  the  case  of  parametric  decay  into  lower  hybrid  wave  at  upper 

hybrid  resonance  with  a linearly  polarized  pump.  (S.  Hiroe  and 

22 

H.  Ikegami's  experiment.  ) 


T = 4 eV 
e 


T.  = 300  K 
1 


N = 3.  4 m torr 


4 - 1 

w = 27r  X 3.  5 X 1 0 sec 
s 


9 - 1 

w.  = 27r  X 4.  1 X 10  sec  = w 
i o 


B = 1350  Gauss  17.  = 27r  x 1 . 05  x 1 o'^sec 17  = 27r  x 3.  85  x 1 O^sec’^ 

O 1 © 
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m.  = 200  X 1.  67  X 10  g 
1 ® 


From  (4.  49),  the  electric  field  associated  with  the  threshold  power  can  be 
calculated  and  the  result  is 


= 2.  25  V/cm 


and  the  experimentally  measured  threshold  field  is  about  2 V/cm. 
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V.  DESCRIPTION  OF  EXPERIMENT  AND  RESULTS 


9-13  20  23  24 

Although  the  theories  and  experiments  ’ ’ of  parametric 

coupling  between  longitudinal  electron  plasma  waves  and  ion  acoustic 
waves  in  an  unmagnetized  plasma  have  been  treated  extensively  by  various 
authors,  only  recently  did  experiments  ’ ’ ’ ^ in  magnetized  plasma  re- 

port on  the  decay  instability.  Porkolab  and  his  colleagues  have  done  a 
series  of  experimental  studies  of  plasma  heating  due  to  the  parametric 
decay  instability  of  plasma  waves  in  a magnetic  field  and  a high-frequency 
electric  field.  These  decay  waves  have  relatively  short  wavelengths  and 
broad  spectra  such  that  quasi-electrostatic  modes  may  be  assumed  and  an 
anomalous  heating  process  can  be  achieved.  In  our  experiment,  we  have  a 
relatively  small  sized  plasma  beam.  Working  on  the  slab  model,  the  ex- 
cited waves  are  of  standing  wave  type  with  well  defined  boundary  condi- 
tions. Therefore,  the  frequency  spectra  of  the  unstable  waves  are  sharp- 
ly spiked.  The  object  of  the  present  work  is  to  study  both  experimentally 
and  theoretically  the  threshold  and  saturation  of  the  parametric  decay  in- 
stability of  the  second  harmonic  of  the  electron  cyclotron  wave  and  ion 
acoustic  wave. 

1.  Experimental  Apparatus  and  Procedure 


The  experiment  is  performed  in  a hollow-cathode-arc-discharge 
25 

(HCD)  plasma  source  , 15  cm  in  diameter  and  a vacuum  chamber  2m 
in  length,  as  shown  in  Figure  1.  The  2m  long  stainless  steel  vacuum 
chamber  is  separated  into  two  sections,  a source  region  and  a drift  re- 
gion, by  a baffle.  In  order  to  get  a highly  ionized  amd  confined  argon 
plasma  beam  in  the  drift  region,  the  magnetic  field  in  the  source  region 
is  a mirror  field,  hence  it  is  hard  to  get  a totally  quiescent  plasma.  In 


FIGURE  1:  SCHEMATIC  DIAGRAM  OF  PINY  HCD  SYSTEM 


our  experiment  a new  procedure  is  introduced  to  create  an  ultrastable 
microwave  sustained  plasma.  An  antenna  horn  is  mounted  on  the  anode 
and  to  send  the  high  power  microwave  (150-200  watts)  into  the  HCD  plas- 
ma beam.  Increasing  the  drift  magnetic  field  and  background  gas  pres- 
sure to  a suitably  high  range  (1.  5 KG-2KG,  1 jx  - 1 • 5fx  of  mercury  pres- 
sure), we  can  then  shut  off  the  source  plasma  completely  by  setting  the 
cathode  gas  flow  and  baffle  current  to  zero.  Yet  the  plasma  still  exists 
in  the  drift  region  sustained  by  the  input  microwave.  This  plasma  is 
very  quiescent  and  well  confined  in  a beam.  Adjusting  the  experimental 
parameters  (such  as  drift  magnetic  field,  drift  gas  pressure,  microwave 
power,  etc.  ) to  proper  values  to  achieve  optimum  operating  conditions, 
we  obtain  a spectrum  as  in  Figure  3,  showing  the  instability  and  its  har- 
monics due  to  parametric  excitation. 

A simplified  block  diagram  of  the  experimental  apparatus  is  shown 
in  Figure  2.  A ffequency  stabilized  Klystron  unit,  operated  in  x-band  is 
used  as  the  microwave  source.  For  the  measurements  of  the  growth  and 
decay  times  of  the  excited  waves,  a fast  PIN  diode  switch  having  20  nano- 
second rise  and  fall  times  is  used  to  modulate  the  microwave  signal  which 
is  fed  into  a Klyst/on  amplifier,  capable  of  1 K watts  cw  output  at  the  fre- 
quency of  9..  23  GHz.  Between  the  PIN  diode  and  the  source  two  variable 
attenuators  are  used  to  adjust  the  output  levels. 

Microwave  power  is  transmitted  through  waveguide  into  the  plasma 
by  means  of  a horn.  Since  the  EM  wave  in  the  waveguide  is  TE^q  mode, 
the  electric  field  in  the  chamber  is  linearly  polarized  in  the  direction 
perpendicular  to  the  uniform  magnetic  field  and  almost  uniformly  distri- 
buted at  the  cross-section  of  the  chamber  near  the  horn.  The  diagnostic 
devices  employed  in  this  experiment  include  axially  and  radially  movable 
Langmuir  probes. 


FIGURE  2:  SIMPLIFIED  BLOCK  DIAGRAM  OF  EXPERIMENTAL  APPARATUS 


FIGURE  3;  FREQUENCY  SPECTRUM  OF  ELECTROSTATIC  ION  CYCLO 
TRON  WAVE 

First  Peak  is  the  Zero  Reference,  Center  Frequency  = w / 2*’ 
= 228  KHz,  P = 1.  2^t,  f = 9.  23  GHz,  P = 173  W (cw).  Line  a 
Vertical  Scale,  90  KHz/Division  Horizontal  Scale 


TIME  DOMAIN  BEHAVIOR  OF  ELECTROSTATIC  ION 
CYCLOTRON  WAVE 

P=  1.  2fi,  f = 9.  23  GHz,  P = 173  W (cw).  Horizontal 
Scale  = 2u8/^m,  Vertical  Scale  = . 5 v/cm 
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In  order  to  parametrically  excite  the  electrostatic  ion  cyclotron 
waves  efficiently,  the  external  dc  magnetic  field  is  suitably  adjusted  such 
that  the  frequency  of  the  pump  field  is  near  or  just  above  the  harmonic  of 
electron  cyclotron  frequency.  We  find  that  in  the  HCD  device,  the  plasma 
is  very  unstable  due  to  the  existence  of  various  types  of  low  frequency  un- 
stable modes  which  are  suspected  to  be;  resistive  drift,  ion  acoustic,  and 
electrostatic  ion  cyclotron  instabilities.  It  is  believed  that  the  energy 
source  responsible  for  the  onset  of  these  instabilities  must  come  from  the 
excess  of  the  free  energy  contained  in  the  plasma  not  at  thermodynamic 
equilibrium.  In  this  regard,  we  find  the  possible  deviations  from  equili- 
brium occurring  in  the  velocity  space  are  due  to  the  mirror  field  in  the 
source  region,  and  not  to  the  configuration  space.  This  conclusion  follows 
from  the  fact  that  all  the  low  frequency  unstable  modes  disappear  after  the 
source  plasma  is  turned  off  completely. 

Strongly  enhanced  signals  at  the  acoustic  wave  frequencies  are  ob- 
served only  when  the  pump  power  exceeds  a threhold  level.  This  is  the 
characteristic  of  parametric  excitation.  With  the  background  gas  pres- 
sure readjusted  to  about  1.2^  the  electrostatic  ion  cyclotron  wave  be- 
comes most  coherent  in  time  domain  as  seen  on  a scope  (Figure  4),  or 
becomes  sharpest  in  the  frequency  domain  as  seen  on  a frequency  spec- 
trum analyzer  (Figure  3). 

2.  Measurement  of  Growth  and  Decay  Times  of  Electrostatic  Ion  Cyclo- 
tron Waves 

The  growth  and  decay  times  of  parametrically  excited  electro- 
static ion  cyclotron  waves  are  measured  as  a function  of  microwave  pump 
power  by  modulating  the  pump  as  illustrated  in  Figure  5.  During  time 
Tj^,  the  power  is  at  level  (threshold  power)  such  that  the 
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FIGURE  5;  GROWTH  AND  DECAY  OF  ELECTROSTATIC  ION  CYCLO- 
• TRON  WAVES  AS  A FUNCTION  OF  PUMP  POWER 
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electrostatic  ion  cyclotron  wave  grows  to  a finite  amplitude  as  shown  in 
the  sketch.  At  the  end  of  time  and  the  beginning  of  T the  power  is 
switched  to  a new  level  P2.  If  P2  < the  wave  decays  to  zero  during 

time  T 2’  providing  that  T2  is  of  sufficient  duration,  as  shown  in  Figure 
5a.  If  P2  ^ ^th  ^ ^2’  wave  decays  during  from  the  am- 

plitude excited  by  P^  to  a new  nonzero  steady  state  amplitude  corres- 
ponding to  P2  as  shown  in  Figure  5b.  Modulating  the  power  as  described, 
the  growth  and  decay  times,  and  the  amplitlide  of  the  electrostatic  ion  cy- 
clotron wave,  are  measured  as  a function  of  P2.  for  both  below  and  above 
the  threshold.  The  total  time  + T2  of  the  microwave  signal  is  appro- 
ximately 2.  3 milliseconds  and  is  typically  7 5%  of  this  period. 

Photographs  containing  many  traces  of  the  electrostatic  ion  cyclo- 
tron waves  are  used  to  obtain  the  growth  and  decay  time.  Since  the  plas- 
ma beam  is  isolated  by  a large  vacuum  chamber,  the  signal  detected  by 
the  optical  system  (through  optical  fiber  bundles  to  the  photomultiplier 
tube)  is  too  weak  to  overcome  the  inherent  noise  of  the  photomultiplier 
tube.  Due  to  this,  a Langmuir  probe  to  detect  the  electrostatic  ion  cyclo- 
tron waves  is  used.  The  results  thus  obtained  are  fairly  good  even  with 
the  slight  difficulty  of  synchronizing  the  excited  signal  with  the  modulating 
signal.  This  probe  detection  technique  is  also  used  to  measure  the  steady 
state  amplitude  of  the  potential  oscillation,  thus  enabling  us  to  obtain  in- 
formation on  the  saturation  mechanism. 

A series  of  photographs  containing  the  average  growth  and  decay 
of  the  electrostatic  ion  cyclotron  wave  is  shown  in  Figure  6 where  the 
pump  power  is  modulated  as  in  Figure  5b  with  T^  = 1.7  ms  and 
T2  = . 6 ms  and  both  P^  and  Pg  are  above  the  threshold  level.  First, 

Pj  is  set  to  a level  such  that  the  instability  is  most  coherent  with  the 


(c)  Pg  = 65.  5 W,  Pj  = 173  W 


FIGURE  6:  DECAY  AND  GROWTH  OF  ELECTROSTATIC  ION  CYCLOTRON 

WAVE  FROM  ONE  STEADY  STATE  AMPLITUDE  INTO  ANOTHER 

P=  1.  2li,  f = 9.  23  GHz,  Horizontal  Scale  = . 2 ms/cm 
o 

To  = • 6 ms,  T,  = 1.7  ms 


FIGURE  7:  PHOTOGRAPH  SHOWS  THAT  THERE  IS  A TRANSITION 
PERIOD  BEFORE  WAVE  REALLY  STARTS  TO  GROW 
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lowest  possible  neutral  gas  pressure.  We  then  progressively  increase 
the  voltage  amplitude  of  the  modulating  signal  (which  is  applied  to  the 
PIN  diode  switch)  to  lower  the  P2  level,  and  the  series  of  photographs 
are  taken  with  the  scope  trace  triggered  to  start  when  the  power  is 
switched  to  P,j  level.  The  initial  decay  rate  of  the  instability,  when 
power  is  switched  from  P^^  to  P2,  is  difficult  to  be  calculated  because 
of  its  strong  dependence  on  the  initial  amplitude  of  the  instability;  there- 
fore, the  average  decay  rate  as  shown  in  Figure  8 is  calculated  by  aver- 
aging over  the  e -folding  time,  i.  e.  , the  time  it  takes  for  the  amplitude 
decreases  to  l/e  of  its  initial  value.  Figure  6 also  shows  the  existence 
of  a transient  before  the  wave  starts  to  grow.  This  may  be  explained  by 
noting  that  during  both  T^  and  T^  of  the  growth  experiments,  the  electron- 
neutral  collision  frequency  is  much  larger  than  the  average  decay  rate  of 
the  electrostatic  ion  cyclotron  wave  shown  in  the  Figure  8.  During  T2 
the  electron  temperature  of  the  plasma  is  lowered  because  of  decrease  in 
pump  power,  and  at  the  end  of  T2,  the  electron  temperature  rises  back 
quickly  to  its  original  level,  but  it  takes  finite  time  for  the  frequency  of 
the  instability  to  switch  back.  A typical  multiple  trace  photograph  shown 
in  Figure  7 clearly  displays  this  kind  of  delay  phenomenon.  Nevertheless, 
the  linear  relation  between  the  initial  growth  rate  y and  is  clearly 
demonstrated  in  Figure  9.  The  curve  in  Figure  8 shows  that  as  P2  pro- 
gressively decreases,  (as  noted  in  the  series  of  photographs)  the  decay 
rate  increases.  Eventually  P2  reaches  a power  level  such  that  it  takes 
the  wave  the  period  T2  to  decay  to  zero  as  shown  in  Figure  10,  and  this 
corresponds  to  the  threshold  power. 

Decreasing  T2  and  increasing  the  total  period  T^  + T2  (i.  e.  , de- 
creasing the  modulation  frequency)  of  the  microwave  signal,  it  is  possible 


4 


FIGURE  9:  INITIAL  GROWTH  RATE  OF  ELECTROSTATIC  ION  CYCLOTRON 
WAVE  FROM  ONE  STEADY  STATE  AMPLITUDE  TO  ANOTHER  AS 
A FUNCTION  OF  POWER  P.,:  P = 173  W,  P = 1 . 2 u,  f = 9.. 23  GHz 


r 


I 
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FIGURE  10;  DECAY  OF  ELECTROSTATIC  ION  CYCLO- 
TRON WAVE  FROM  STEADY  STATE  AM- 
PLITUDE INTO  NOISE  FOR  PUMP  POWER 
JUST  BELOW  THRESHOLD  LEVEL 

Pg  = ^7.  1 W.  P^  = 141  W,  Tg  = . 2 ms, 

Tj  = 2.  1 ms,  Horizontal  Scale  = . 1 ms/cm 
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to  lower  the  power  level  P2  to  almost  zero.  Single  and  multiple  traces 
photographs  shown  in  Figure  11  display  the  growth  (and  decay)  of  the  elec- 
trostatic ion  cyclotron  wave  as  power  P2  is  varied  from  0 watt  to  173 
watts.  Only  a portion  of  the  total  period  is  shown  in  the  photograph.  The 
wave  starts  to  grow  only  with  the  pump  power  exceeding  the  threshold 
level,  and  it  becomes  an  eigenmode  of  the  collisional  plasma.  During 
the  initial  buildup  of  the  potential  oscillation,  the  energy  density  arising 
from  spontaneous  emission  of  the  background  plasma  is  comparable  to  the 
energy  density  produced  by  the  parametric  interaction.  Therefore,  to 
determine  the  initial  growth  rate  we  must  calculate  the  growth  time  only 
after  the  wave  has  a small  finite  amplitude.  This  value  is  also  shown  in 
Figure  9,  and  is  consistent  with  the  theoretical  calculation  from  equation 
(4.  17). 

There  are  two  kinds  of  damping  in  our  system.  Linear  or  natural 
damping  occurs  when  the  pump  is  below  the  threshold  level.  Once  the 
threshold  level  is  exceeded,  the  amplitude  of  the  instability  would  theo- 
retically grow  to  infinity.  However,  due  to  the  nonlinear  damping,  the 
wave  always  reaches  a saturated  state.  Therefore,  the  nonlinear  damp- 
ing would  exactly  balance  the  initial  growth  rate  of  the  wave  at  that  parti- 
cular power  level.  If  the  pump  power  is  suddenly  reduced  to  zero  as  in 
the  case  in  Figure  lib,  its  initial  decay  rate  must  be  equal  to  the  sum  of 
the  linear  damping  rate  and  the  initial  growth  rate  of  Figure  11a,  Hence 
we  have  a linear  damping  rate  of  about  0.  6 x lO^sec  ^ for  the  electro- 
static ion  cyclotron  waves;  this  value  is  consistent  with  the  value  calcu- 
lated from  the  tail  of  the  decay  diagram.  From  this  value,  we  can  con- 
clude that  the  linear  damping  mechanism  is  due  to  ion-neutral  crllisions 
with  the  ion  temperature  of  about  0.  6 eV.  This  value  of  ion  temperature 


(a)  SINGLE-SHOT  PHOTOGRAPH 
Pg  = 0 W,  Pj  = 173  W,  P = 1.  2 
Horizontal  Scale  = 50  ns/cm 


1 

(b)  MULTIPLE -TRACE  PHOTOGRAPH 

Pg  = 0 W,  Pj  = 173  W,  P=  1.  2 fi 

Horizontal  Scale  = . 1 ms/cm 

1 



has  been  verified  experimentally  in  Kristal's  ” thesis  by  using  a Fabry- 
Perot  interferrometer. 


3-  Frequency  of  Electrostatic  Ion  Cyclotron  Wave  as  a Function  of  Pump 
Power 

The  frequency  shift  of  the  electrostatic  ion  cyclotron  wave  due  to 
pump  power  has  also  been  investigated.  These  experiments  are  con- 
ducted with  cw  microwave  excitation,  and  the  frequency  is  measured  on 
the  low  frequency  spectrum  analyzer. 

Electrostatic  ion  cyclotron  frequency  is  shown  as  a function  of 
power  in  Figure  1 2 for  a fixed  pressure.  As  shown  in  Table  1,  the  fre- 
quency shift  is  linearly  proportional  to  the  change  of  electron  temperature 
with  some  discrepancy.  This  can  be  explained  from  the  linear  parametric 
theory,  which  indicates  a frequency  shift  due  to  finite  growth  rate.  After 
substituting  all  the  known  parameters  into  equation  (4.  20),  this  discrep- 
ancy becomes  more  evident.  However,  with  a minor  correction,  the 
linear  relationship  between  the  frequency  shift  and  electron  temperature 
change  is  clearly  shown  in  Figure  13. 

Electron  density  with  the  presence  of  the  pump  wave  is  also  mea- 
sured by  using  the  probe  technique.  As  shown  in  Figure  14  the  electron 
density  at  the  center  of  the  plasma  column  decreases  as  the  instability 
increases.  But  the  pump  power  increases  over  some  level,  the  ionization 
rate  may  cover  the  enhanced  diffusion  rate,  causing  the  density  to  increase 
again. 


4.  Summary  of  Experimental  Parameters 


Table  2 summarizes  the  basic  experimental  parameters.  The  values 
given  are  for  a magnetic  field  of  1.  64  kilograms  and  for  a pressure  of  1.2^ 


TABLE  1 --  EXPERIMENTAL  AND  THEORETICAL  DATA  IN  SEVERAL  PUMP  POWER  LEVELS 


density 


Af,  (X 


ATg(eV) 

FIGURE  13:  MEASURED  ELECTROSTATIC  ION  CYCLOTRON  FRE- 
. QUENCY  SHIFT  AS  A FUNCTION  OF  ELECTRON  TEM 
PERATURE  CHANGE 


FIGURE  14;  PLASMA  DENSITY  AS  A FUNCTION  OF  PUMP  POWER 


-59- 


1 


TABLE  2 - EXPERIMENTAL  PARAMETERS 


Quantity 

Symbol 

Value 

Ion  (Electron)  Density 

n 

8.8x10^^/  cm^ 

Neutral 

KJ 

N 

3.9x10^^/ cm^ 

Ion  Temperature 

T. 

1 

0,  6 eV 

Electron  Temperature 

T 

6 eV 

eBo 

e 

Ion  Cyclotron  frequency 

1 

3. 95  X 1 0^/ sec 

eB 

2.  9 X 1 0^  *^/sec 

Electron  cyclotron  frequency 

e 

Electrostatic  Ion  cyclotron  frequency 

P k^T  /M 

0) 

s 

1.3x1  0^/ sec 

Second  harmonic  of  electron  cyclo- 
tron frequency 

" i 

5.  8 X 1 0^  *^/sec 

Pump  frequency 

CO 

0 

5.  8 X 10^^/sec 

Ion-neutral  collision  frequency 

u,  =2r 

\n  s 

1.  2 X 10^ /sec 

Electron-neutral  collision  frequency 

^en  ~ 

5.2x10  /sec 

. 2 l/2 

4Tn  e 

1 . 95  X 1 0^/ sec 

Ion  plasma  frequency  ( — ) 

u 

pi 

. 2 l/2 

47rn  e 

5.3  X 1 0^  *^/sec 

Electron  plasma  frequency  ( — ) 

U) 

Pe 

T 1/2 

Debye  length  ( „ ) 

^D 

1 . 94  X 1 0 ^ cm 

4Tn  e 
o 

Radial  scale  length 

X 

r 

. 4 8 cm 

Wave  number  ^ 

k 

3.  3/ cm 

r 

8T.  1/2 

Ion  Larmor  radius  ( /n. 

ttM  1 

"^Li 

. 48  cm 

00 

H 

to 

Electron  Larmor  radius  ( ) / n 

r_ 

.57x10  cm 

irm  e 

Le 

3T.  l/2 

Ion  thermal  speed  (“j^  ) 

^Ti 

2.  1 X 1 0 cm/ sec 

3Tg  1/2 

Q 

1.8x10  cm/  sec 
- 23 

6.  68  X 10  gram 

Electron  thermal  speed  ( ) 

^ m 

Ion  mass 

^Te 

M 

Electron  mass 

m 

0.  91  X 1 0 gram 

L 
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Note  that  V.  includes  charge  transfer,  which  accounts  for  about 
in  “ 

- 1 4 2 

half  of  the  total  collision  rate.  A cross-section  of  1 . 3 x I 0 cm  is 
used  in  the  calculations. 


5.  Evidence  of  Parametric  Excitation  of  Electrostatic  Ion  Cyclotron 
Wave 

Experimental  evidence  and  data  are  presented  here  to  show  that  the 

electrostatic  ion  cyclotron  waves  are  indeed  parametrically  excited  by 

the  microwave  pump.  It  is  found  that: 

a.  The  polarization  of  the  excited  low  frequency  wave  is  in  the  direc- 
tion perpendicular  to  the  dc  magnetic  field. 

b.  It  is  a standing  wave. 

c.  The  excited  low  frequency  wave  is  the  electrostatic  ion  cyclotron 
wave. 

d.  The  excited  waves  start  to  grow  only  when  the  microwave  pump 
reaches  threshold  power. 

e.  Its  growth  and  decay  rates  follow  the  theoretical  linear  dependence 
with  power  as  predicted  by  parametric  theories. 

a.  First  the  axially  movable  probe  is  calibrated  with  respect  to  the 
radially  movable  one.  Then  both  probes  are  used  to  measure  the 
field  strength  of  the  wave  at  the  same  point  but  in  two  different 
polarizations.  We  find  that  the  axial  field  is  much  smaller  than 
the  field  component  in  the  direction  transverse  to  the  dc  magentic 
field. 

b.  Using  the  radially  movable  probe  to  measure  the  amplitude  of  the 
transverse  field  as  a function  of  the  radial  position,  the  amplitude 
is  found  to  be  maximum  at  the  center  and  reduces  to  zero  just  out- 
side the  plasma  beam.  Hence  we  conclude  this  excited  mode  is  of 
standing  wave  type,  and  its  wave  number  may  be  calculated  from 
the  diameter  of  the  plasma  beam. 

c.  To  support  the  contention  that  the  excited  low  frequency  wave  is 
electrostatic  ion  cyclotron  wave,  we  note  that  the  observed  fre- 
quency is  almost  independent  of  the  plasma  density  and  dc  magnetic 
field;  therefore  the  possibility  of  exciting  lower  hybrid  wave  or  har- 
monics of  ion  cyclotron  wave  is  excluded.  We  also  note  that  the 
observed  electrostatic  ion  cyclotron  frequency  matches  the  one  de- 
termined by  the  linear  dispersion  relation  with  known  wavenumber, 
electron  temperature,  ion  mass  and  dc  magnetic  field.  Since  the 
linear  dispersion  relation  of  the  electrostatic  ion  cyclotron  wave  is 
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J = nf  +k^  jfji 

the  frequency  shift  due  to  the  change  of  electron  temperature  may 
be  expressed  as 

1 

This  linear  relationship  has  been  shown  in  Figure  13. 

d.  Multiple  trace  photographs  shown  in  Figure  6 and  Figure  10  display 
both  decay  and  growth  phenomena  of  the  electrostatic  ion  cyclotron 
wave.  Decay  waves  have  finite  steady  state  amplitude  only  for  pump 
powers  above  a threshold. 

e.  Possibly  the  most  convincing  evidence  of  parametric  excitation  of  the 
electrostatic  ion  cyclotron  wave  lies  in  the  main  theme  of  this  effort, 
i.  e.  , the  growth  and  decay  rates  of  the  excited  ion  acoustic  waves 
are  linearly  related  to  pump  power  as  shown  in  Figures  8 and  9. 

6.  Simplified  Theoretical  Explanation 


The  linear  parametric  theory  derived  in  the  preceding  chapter  gives 
us  a consistent  prediction  of  threshold  power  and  initial  growth  rate. 

Since  it  does  not  include  the  saturation  effects,  it  fails  to  describe  the 
growth  or  decay  of  the  unstable  wave  from  one  saturated  amplitude  to 
another.  Nevertheless,  the  result  given  in  Eq.  (4.  19)  leads  us  to  pro- 
pose a phenomenological  model  for  wave  growth  and  decay  in  the  presence 
of  saturation  effects  as  the  pump  power  is  near  threshold.  The  amplitude 
A of  the  electrostatic  ion  cyclotron  wave  may  be  governed  phenomeno- 
logically as 


A - C3  A 


n+1 


(5.  1) 


where  is  the  square  of  the  electric  field  in  the  plasma  produced  by 

2 

the  pump  and  E is  the  square  of  the  threshold  field.  F is  the  linear 
c s 
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damping  rate  and  is  the  phenomenological  nonlinear  damping  rate 

and  n is  an  integer  of  unknown  value  which  may  be  determined  from  the 
experimental  data. 

The  exact  solution  of  (5.  1)  for  n > 0 is 


A(t)  = 


A e 
o 


yt 


I - 1) 

y 


l/n 


(5.  2) 


where  y = F 


- 1 


is  the  linear  growth  (or  decay)  rate  and  A^=  A(0) 


is  the  initial  amplitude  of  the  wave. 
As  t -*  00  „ 


- J£ c 


>0  for  y > 0 


(5.3) 


and 


A(oo)  = 0 


for  y < 0 


(5.  4) 


Hence  the  value  of  n can  be  experimentally  determined  by  plotting  the  n^^ 
power  of  the  steady  state  amplitude  of  the  electrostatic  ion  cyclotron  wave 
as  a function  of  pump  power  for  y > 0 as  shown  in  Figure  15.  With  n = 2, 
this  graph  becomes  linear,  and  is  in  agreement  with  the  prediction  of  (5.  3). 
Setting  n = 2 in  (5.  2)  and  (5,  3),  we  obtain 


A{t)  = 


A e 
o 


yt 


1 1—2^  (gSyt  _ j 


1/2 


(5.  5) 


and 


A(«.)  = (f-) 

3 


1/2 
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1/2 


for  y > 0 


(5.  6) 


respectively. 

Now  we  proceed  to  explain  the  experimental  results  shown  in  Fig- 
ures 8 and  9 by  using  (5.  1). 

Assuming  the  pump  power  is  switched  from  to  P2  (or  P2  to  P^ ) 
as  in  the  case  shown  in  Figure  5b,  then  the  wave  decays  (or  grows). from 
original  steady  state  amplitude  A(0)  to  another  steady  state  amplitude  A(°o)- 

Rewriting  Eq.  (5.  1 ) as 


dA 

dt 


e2 

E 1 2 

£ 

A = 

Tr2  ^2. 

^1  ' ^2  2 2 
rs(  ) + c3(a2-a2) 

c 

C 

(5.7) 


for  the  growth  case 


where  and  E2  are  the  field  corresponding  to  P^^  and  P2  respec- 

Eg 

r (-1  - 1) 

2 " ^c 

tively,  and  A^  = is  the  steady  state  amplitude  of  electro- 

3 

static  ion  cyclotron  wave  at  P^  power  level,  and 


^2 

2 2 1 

dA  _ 
dt 

^2  2 

A=  - 

^1-^2  2 2 

c 

c 

decay  case 


A for 


(5.  8) 


! *%  Mi** 


ion  cyclotron  wave  at  power  level. 

Then  Eqs.  (5.  7)  and  (5.  8)  show  that  the  initial  growth  rate  (or  de- 
cay rate)  is  linearly  proportional  to  the  power  level  P2‘ 
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VI.  NONLINEAR  SATURATION  MECHANISM 


1.  Introduction 

We  have  so  far  considered  only  linear  paranaetric  theory  and  a phe- 
nomenological nonlinear  equation  which  describes  the  low  frequency  unstable 
wave.  Since  saturation  phenomenon  is  always  present  in  experiments,  an 
important  question  must  be  raised  in  order  to  describe  the  kind  of  nonlinear 
damping  mechanisms  which  cause  the  instabilities  to  saturate.  There  have 

been  earlier  attempts  to  describe  the  mechanism  as  due  to  induced  scatter- 
27  28  29  30  31 

ing  ’ , and  cascading  ’ ’ of  the  side  band  into  even  lower  frequency 

32 

waves  and  resonance  broadening  . Other  possible  mechanisms  include 
pump  depletion  and  quasi-linear  effects  (i.  e.  , heating).  There  is  no 
generalization  to  be  applied,  because  the  nonlinear  effects  may  vary  for 
different  experimental  systems.  As  an  example,  if  one  nonlinear  damping 
mechanism,  among  other  competing  phenomena,  becomes  dominant  and 
cause  to  saturate  before  the  pump  depletion,  then  we  must  rule  out  pump 
depletion  due  to  priority,  or  if  the  unstable  wave  is  coherent,  then  we  must 
reject  the  mechanism  due  to  the  cascading  of  the  sideband.  Therefore, 
the  dominance  among  different  mechanisms  depends  greatly  on  the  geome- 
try of  the  individual  system,  and  the  mode-types  involved.  Disregarding 
the  effect  of  inhomogeneity,  systems  may  be  classified  into  three  classes: 

(1)  Infinite  plasma 

When  the  dimensions  of  the  plasma  is  much  larger  than  the  wave- 
length of  the  instabilities,  boundary  conditions  need  not  be  considered. 
Therefore,  the  excited  modes  are  traveling  waves,  and  the  spectrum  may- 
be strongly  broaded.  In  such  a case,  the  saturation  of  instabilities  may  be 
due  to  the  pump  being  depleted  to  the  threshold  level.  Another  possibility 
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is  the  interaction  of  the  excited  large-amplitude  oscillations  among  them- 
selves, causing  a continuous  flow  of  energy  toward  the  large  wave  number 
region.  This  oscillation  energy  is  eventually  dissipated  into  heat  by  col- 
lisional  damping,  or  via  collective  effects.  The  threshold  level  of  the  in- 
stability is  thus  increased,  and  a nonlinear  saturated  state  may  eventually 
be  reached. 

(2)  Bounded  plasma  with  sheath 

If  the  spectra  of  the  excited  modes  are  sharply  defined  because  of  the 
confines  of  the  boundary,  the  potential  oscillations  of  the  low  frequency 
waves  may  carry  the  ion  bunch  into  the  sheath  region.  Therefore,  recom- 
bination of  ion  bunches  at  the 'sheath  may  be  the  dominant  nonlinear  damping 
mechanism  in  this  case.  Since  the  recombination  rate  is  linearly  propor- 
tional to  the  density  of  the  ion  bunch,  this  nonlinear  damping  process  cor- 

24 

responds  to  the  n = l case. 

(3)  Bounded  plasma  without  sheath 

Excited  modes  are  of  the  standing  wave  type,  producing  a well  defined 
frequency  spectra.  From  supercritical  stability  theory,  anomalous  diffusion 
of  the  plasma  due  to  the  low  frequency  unstable  oscillations  may  be  one  of 
the  nonlinear  damping  mechanisms,  which  corresponds  to  n=2  case. 

Another  possible  mechanism  is  the  nonlinear  harmonic  (or  subharmonic) 
generation,  with  second  harmonic  generation  corresponding  to  n=  2 case. 

2.  Theory 

In  this  section  the  description  of  the  nonlinear  damping  mechanisms 
for  case  (3),  bounded  plasma  without  sheath,  will  be  given.  The  experi- 
mental results  in  Chapter  V show  that  (a)  a density  decrease  in  the  central 
portion  of  the  plasma  column  coincides  with  an  increase  in  amplitude  of 
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the  instability.  When  pump  power  is  increased  over  a certain  level,  the 
ionization  rate  would  cover  the  enhanced  diffusion  rate,  causing  the  density 
to  increase,  (b)  Harmonics  of  the  instability  are  also  excited.  The  ampli- 
tude relative  to  that  of  the  fundamental  is  small  if  the  pump  level  is  low, 
but  increases  with  an  increase  in  pump  power. 

It  should  be  pointed  out  that  since  the  excited  mode  is  an  electrostatic 
ion  cyclotron  wave,  these  harmonics  are  still  in  the  resonance  region,  and 
they  have  the  same  linear  damping  rate  as  that  of  the  fundamental  mode. 
Neglecting  the  nonlinear  damping  of  the  harmonics,  the  nonlinear  damping 
rate  of  the  fundamental  instability  due  to  its  harmonic  generation  may  be 
defined  as  follows 


2T. 


N 


s 1 


2 

2 

6E  (k  ,w  ) 
s — s 

Aw  , = 2r  { 
si  s 

6E  (2k,  20)  ) 
s — s 

Aw  + 
s 2 

6E  (3k  , 3w  ) 
s — s 

Aw  , + . . } 

s3 


(6.  1) 


or 


Aw 

s2 

Aw  , 
s 1 


6E  (2k,  2w  ) 
s — s 

2 Aw 

+ 

6E^(3k.3w^) 

6E  (k,w  ) 
s'—’  s 

Aw 

s 1 

6E  (k,  w ) 
s — s 

(6.  2) 


where  the  left  hand  side  of  (6.  1)  is  the  effective  instability  energy  deple- 
tion rate,  and  the  right  hand  side  is  the  damping  rate  of  the  energy  of  the 
excited  harmonies  in  steady  state,  (i=l»  2,  3 ...  ) is  the  bandwidth  of 

each  spectrum.  Experimental  observations  indicate  that  such  an  infinite 
series  may  be  truncated  after  the  third  harmonic  when  pump  power  is  not 
too  strong. 

Knowing  that  the  nonlinear  damping  rate  in  steady  state  must  equal 

N 

to  the  initial  growth  rate.  ^ can  be  calculated  from  (6.  2)  and  Figure  3, 


however,  this  value  is  much  smaller  than  the  initial  growth  rate  calculated 
from  (4.  17)  . We  may  conclude  that  there  must  exist  another  kind  of  non- 
linear damping  mechanism  which  dominates  in  the  saturation  process. 

It  is  well  established  that  diffusion  exists  in  a magnetically  confined 
beam  plasma,  and  usually  pure  coherent  longitudinal  modes  does  not  cause 
an  enhanced  plasma  loss  across  the  magnetic  field.  With  the  plasma  dimen- 
sion smaller  than  the  wavelength  of  the  oscillation  in  question,  magnetic 
field  causes  the  charged  particles  to  gyrate  and  move  towards  the  side 
boundary  to  produce  a field  which  causes  the  oscillating  mode  to  become 
elliptically  polarized.  This  induced  field  may  cause  an  enhanced  diffusion. 
While  ordinary  diffusion  of  electrons  and  ions  across  a magnetic  field  is 
caused  by  collisions  with  non- identical  particles,  it  has  been  observed  that 
certain  oscillations  seem  to  enable  the  plasma  to  acquire  an  enhanced  dif- 
fusion rate.  The  behavior  of  a low  temperature  plasma  is  then  largely 
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governed  by  the  diffusion  processes  involved,  but  the  anomalous  diffusion 

can  be  visualized  as  a consequence  of  additional  particle  collision  with  the 

electric  field  of  the  oscillating  instability.  This  can  also  be  regarded  as 

made  of  ^ x ^ drift  motion  of  the  particles,  and  most  readily  understood 

from  the  test  - particle  point  of  view.  Let's  write  the  equation  of  motion  for 

a single  ion  moving  in  a wave  field  with  the  presence  of  a uniform  magnetic 

field  B z as  follows: 
o 


dt  ■ - 


(6.  3) 


1 

® ^ = e^(£(t),  t)+Mf2j^^(r_(t),  t)  X z 


and  the  solution  of  (6.  3)  may  be  expressed  as 
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v(r(t).t)  =R.(t)  . v(0)+-^  / dt'  R.  (t-f).  E(r  (t').t') 

o ~ 

, t 

r(t)  = r(0)  + -^  L.  (t).  v(0)+^  f dt'  L.  (t-f).  E (r(t').  t') 


(6.4) 


where 


i_  A 

n.  dt 
1 


(t)  = 

cos  il.t 

1 

-sinil.t 

sinfl.t 

I 

cos  n.t 

* 

0 

0 

L.  (t)  = 

sinJl.t 

1 

- (1  -cos  n.t ) 

1-cosil.t  0 

sinfl.t  0 

1 

0 

1 

0 

1 

= 1 

1 

0 

1 

0 n.t 

1 

Thus  the  position  R (t)  of  the  guiding  center  is  obtained  as 


R(t)  = r(t)  -j|-zxv(t) 


t t 

= /dt'^(r(t')+z{v  (0)t+^  /dt"E  (r^t").t")}  (6.5) 


and  the  transverse  drift  velocity  of  the  guiding  center  is  given  by 
dR^  (t) 

V^{r(t).t)  = £xE  (r(t),t) 


(6.  6) 


where 


t 

R (t)  = R - /dt'E(r(f),t) 

o o 


(6.7) 


with  the  aid  of  Eq.  (6.  5),  V^(£(t),t)  nniay  be  expanded  with  respect  to  the 
center  position  ^(t)  as 


1 - 


y^(r(t),t)=V^(R(t),t)  + zxv(t)  . Vy^(R(t),t)+  .. 


(6.8) 


= - (R(6)|  t)  + j^zxv  (t).  VE(R(t),t)+.  . . } 

i 


This  would  give  the  following  relation  immediately 


w 
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V (R(t),t)  = - zxE  (R(t),t) 

X B 


where  (t)  may  be  approximately  expressed  as 


-i  /dt'E(R(t'),t') 

o o 

therefore,  in  the  Euler's  coordinate  system  we  have 


(6.9) 


(6.  10) 


V (R.t)  = - ixE  (R,t) 
X — ts  — — 


and 


V^(R.t)  = ^ E (R.t) 


(6.  11) 


(6.  12) 


Due  to  the  spatially  dependent  force  field,  there  is  a force  difference 
between  two  different  points,  i.e.  , finite  excursion  of  the  oscillating  ion 
brings  itself  into  regions  of  different  field  intensity.  The  explicit  expres- 
sion for  the  statement  above  is  given  by 


(6.  13) 


where  6^  = ■ ^2 


now  let's  define  a force  due  to  the  force  difference 

(2) 

W \ (Rr'>  X V5.r‘)>=<5R.v^v^(R,,t)> 

(6.14) 

where  < > is  meant  to  ensemble  average 

t t t' 

5R=-  zx/dt'E(Rj,t')+S{v  (0)t  + ^ /dt'/dt"E  (Rj,t")}  (6.15) 


O O ■■  o o 

and  assume 

E (R  j , t)  = A cos  k . sin  ( w t + <()  ) 


(6.  16) 
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where  k = kx  and  4)  is  a random  phase  which  is  a constant  during  the 
period  of  two  successive  collisions. 

Substituting  (6.  15)  and  (6.  16)  in  eq.  (6.  14),  yields 


(2) 

'x(^l 


dV  (R,,t)  _2  t 


= < - ^^  k A [ J^dt' cos  (k.  Rj)sin(wt'+4) )]  sin(k.^j)cos(oJt +4)  )u> 

o 

c2  2 

= < — ^kA  sink.  RjCos  k.  R j[  cos  (cot+4>  )-cos4>  ]cos(ut+4>  )> 

B ^ 

o 

■? 

.2 


= = kA  sin  2k. R.  (1-cos  ut) 

4b2  y 

o 


(6. 17) 


(2) 

therefore,  the  maximum  contribution  to  (R^)  between  two  successive 

collisions  is 

(2)  .2 


V (R.)  = kA^  sin(2k.R,)  t 

^ ^ 4B^  y ^ 

o 

and  with  the  assumption 
u > > — 


(6. 18) 


(6.  19) 


where  t is  the  collision  time, 
c 


(2) 


Because  of  the  spatial  dependence  of  (^j)»  it  will  give  an  additional 

damping  to  V^(Rj»t)  . 

In  fluid  limit  , Rj  is  not  the  initial  position  of  the  particle  trajectory,  it  is 
an  independent  variable.  Thus 


aRj 

at 


= 0 


(6.  20) 


since 


S c d 

4v^(R,,.)+V^(R,.t).Vt  V (R  tHV  (R  )]=^ 

o ^ 

(6.  21) 


■HHUC. 
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This  may  be  written  as 


o ■'  2B  ' 


cos(2k.R,)T  V (R,,t) 

— — 1 C X — 1 


■=|-dT-^<Sr*l-V£  V£r‘> 

o ^ 


(6.  22) 


where 


^in  “ T, 


V = V.  + cos  (2k.  R,)  t 

eff  m __  2 y 1 c 


(6.  23) 


Hence  at  center,  where  Rj  = 0 and 


u = u.  + k^  / u. 
eff  in  _„2  y in 


(6.  24) 


Assuming  A = ^ A,  where  A is  the  amplitude  of  E and  0 < 1 , 

y ^ 

N 

the  nonlinear  damping  rate  of  the  electrostatic  ion  cyclotron  wave  due 

to  the  anomalous  diffusion  becomes 


k^e^A^/zr^ 

o 


(6.  25) 


3.  Comparison  between  Theory  and  Experimental  Results 

As  mentioned  in  Chapter  V,  the  nonlinear  damping  rate  at  steady  state 
must  equal  the  initial  growth  rate.  From  Figure  9 the  initial  growth  rate 

5 

y = 1.  5x10  /sec  for  = 173  watts  is  obtained.  The  amplitude  of  the 
potential  oscillation  for  Pj  =173  watts  can  b»»  found  in  Figure  4.  Since  the 
amplitude  of  the  potential  oscillation  is  the  integration  of  the  amplitude  of 
the  field  oscillation,  then  A =1.  05xk=3.  5 volts/cm  is  found. 
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First  let's  calculate  the  nonlinear  damping  rate  due  to  the  harmonic 
generation.  From  Figure  3 and  Eq.  (6.  2),  we  find 


32x10  /sec 


(6.  26) 


with  the  aid  of  Table  2 and  Eq.  (6.  25),  we  have 


= 10  e^xlO^/sec 
s2  ^ 


since 


N N 5 

r , + r , = y = 1.5x10'’ 

si  s2 


we  find 


n. 


I = .34  « 


(6.  27) 
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VU,  CONCLUSIONS 

The  physical  origin  of  mode  coupling  mechanism  in  an  infinite  plasma 
is  described  as  electrostrictive  effect.  Based  on  this  physical  viewpoint, 
a general  approach  to  analyze  the  parametric  decay  processes  of  plasma 
waves  is  presented.  A set  of  coupled  mode  equations  is  derived  by  using 
Hamiltonian  approach,  and  the  coupling  coefficients  are  derived  from  the 
collisionless  Boltzman- Vlasov  equation.  The  threshold  for  parametric 
amplification  of  an  electrostatic  ion  cyclotron  mode  (longitudinal  wa-ve) 
and  a harmonic  of  electron  cyclotron  mode  (hybrid  wave)  in  the  presence  of 

external  radiation  of  an  appropriate  frequency  has  been  found. 

34 

An  experiment  is  performed  in  a microwave  sustained  plasma.  The 
properties  of  the  process  are  in  qualitative  agreement  with  the  theory. 
Saturation  of  the  instability  is  also  observed  in  the  experiment,  and  the 
nonlinear  mechanisms  are  found  to  be  anomalous  diffusion  and  harmonic 
generation,  again,  in  agreement  with  the  theoretical  results. 
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APPENDIX  A:  NORMALIZATION  OF  THE  ACOUSTIC  MODES 

The  first  step  in  the  quantization  procedure  is  to  write  down  the 

Hamiltonian  density.  Let's  define  a lattice  for  a homogenous  medium  of 

mass  density  p = n^M  , where  n^  and  M are  respectively  the  ion  density 

and  mass,  then  attach  a dressed  vibrating  ion  to  each  lattice  point.  Let  the 

deviation  of  the  mth  ion's  position  from  its  equilibrium  position  be 

35 

given  by  q^  , we  are  led  to  write  for  the  acoustic  modes 


H =1  X 

m=l 


+ Ct  (q  , , -q  ) 

L ^m+1  ^m 


(A.  1) 


Where  for  the  rn=n^  term,  q^  .(.j  = because  of  the  periodic  boundary 

° 1 "3 

conditions  and  P = Mq  , x = m(  — ) and  C.  is  the  force  constant, 

m m m n L 

o 

The  first  term  in  the  brackets  corresponds  to  the  kinetic  energy  of  the 


oscillating  ions,  the  second  to  the  potential  energy  associated  with  the 
Hooke's  law  force.  We  shall  now  proceed  to  rewrite  H so  that  it  has  the 
same  form  as  the  Hamiltonian  for  the  linear  harmonic  oscillator.  Let  us 
express  the  coordinates  q^  and  momenta  P^  in  terms  of  the  traveling 
wave  normal  mode  expansion 


(k)  ,1  .2  , S-m,  , ^ (k) 

<1  = ( — xTTTT  ) (a-i  e ) and  q = Z <3 

^m  2n  Mw ' ' ' k ^m  , ^m 
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(A.  2) 
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m ^m  ' 2n  ' k k 
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where  we  have  a,  = 
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'k'®  ^ 2 ^^k^k*  ^kH  ^k  \ ^k  ^ 
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< V+ 1 - ■’m'  ' ^ TC-  <=‘k\  + \>k  S “k  + ^k  *k  ' 

A»1  K Xj 


. • . H = j Z “k  ^^k^k  ■'■  ^k  ^k  ^ 
k 


(A.  4) 


(A.  5) 


(A.  6) 


If  and  a^  are  now  defined  in  terms  of  the  new  variables  and 

as  follows: 


, M ,2  - ,o(k)  . 


a+  = (M  )2(^.  Q(k).i  K^i) 
k ' 2w'/  ' k ^ M ' 

k 


then  above  Hamiltonian  reduces  to  the  harmonic  oscillator  form 


(A.  7) 


(A_8) 


„ 1 ^ r . ,,  ,2^(k)^n 


(A.  9) 


With  the  Hamiltonian  in  the  form  above  we  proceed  to  quantize  the  acoustic 

(k)  (k) 

vibrations  by  postulating  that  the  Q'  , K are  operators  that  satisfy  the 
poisson  bracket  relations  { ^ } = ^kk'  ' 

Now  let  us  normalize  the  field  produced  by  the  acoustic  vibrations 
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=-E  E 
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.-.  5E^^^=-  ^ 


M 1 2 Q(k) 

e In  k 


(A.  12) 


Where  is  the  resonant  frequency  of  the  medium  without  pump  field. 

If  the  pump  field  presents  in  the  medium  and  near  threshold  level  or  above, 

the  resonant  frequency  will  shift  to  undamped  frequency  w , in  this  case 

k 

we  can  keep  the  unperturbed  Hamiltonian  H as  before,  but  change  «'  to 

o ^ k 

for  the  relation  between  6E^^^  and  . 
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It  is  easy  to  be  generalized  to  three  - dimensional  case 


„(k)  _ 1 y 

^o  - 2 


— ^ + M w'2 

M k cr 


(A.  15) 


6e(^U.M.  JL_  ^2  (k) 


(A.  16) 


We  can  also  generalize  above  procedure  to  an  anisotropic  case,  for 

example,  if  there  is  a uniform  magnetic  field  B in  the  z direction 
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.2 


(k) 


(A.  17) 
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— irT e e fn 

y o y o 
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eB 


where  we  have  { =6  6 and  Q.  = -r^ 

a 8 kk'  ajS  i Me 
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APPENDDC  B:  NORMALIZATION  OF  THE  OPTICAL  MODES 


Because  of  the  big  difference  between  the  masses  of  electron  and  ion, 
ion  cannot  follow  the  high  frequency  electron  oscillation,  hence  we  can  con- 
sider this  kind  oscillations  are  optical  modes  in  which  electrons  and  ions 


vibrate  against  one  another.  Let's  expand  the  displacements  of  the  individ- 


k 

ual  electrons  in  terms  of  the  normal  coordinates  U and  A,  , A, 

k k 
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In  terms  of  these  normal  coordinates  we  may  write  the  total  energy  at  any 
time 
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Now  that  we  have  obtained  the  Hamiltonian  classically  we  may  move  directly 
to  canonical  operator  form.  Then»^^^^  and  become  canonical  variables 
satisfy  the  Poissan  bracket  relation 

The  final  stage  is  to  normalize  the  field  produced  by  the  optical  vibrations 
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For  the  same  reason,  if  pump  field  is  present  and  near  or  above  the  thresh- 
old level 
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In  three  dimensional  case  we  have 
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We  can  also  generalized  above  procedure  to  the  case  that  there  is  an 
uniform  magnetic  field  in  the  z direction 
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jJRITy  cC  AS3I  “ 1 C a T t C S Oa'a  Enivf’f-i/ 


'^The  plasma  .3  pumped  by  a m.oaochrom.atic  electromagnetic  wave  with 

frequency  which  decays  into  t"vo  decay  waves  w'^'^  and  ^ (e.  g.  , a 

plasm.on  and  a phonon) .with  cj  =uj  . Either  one  of  the  -olasm.a  modes 

y o 1 S f[ 

coupled  with  the  pum.p  wave  (through  the  linearized  Vlasov  ecuaiion)  to 
induce  a polarization  which  m*ay  act  as  the  source  of  another  m.ode.  A 
general  procedure  to  calculate  the  induced  polarizations  is  introduced  by 
transformiing  the  linearized  Vlasov  equation  into  the  oscillating  frame  of 
reference,  in  which  the  equation  can  then  be  solved  by  the  method  of 
characteristics.  The  threshold  power,  initial  gro%v;h  rate  and  the  frequency 
shift  (due  to  finite  initial  growth  rate)  of  the  paramietric  instabilities  can  be  ' 
obtained  from,  the  coupled  miode  equations  and  the  results  obtai.ned  can  be 
i applied  to  all  the  wavele.ngth  regio.n.  Several  interesti.ng  cases  are  pre- 
sented and  com.pared  favorably  with  experim.ents. 

The  simultaneous  excitation  of  second  harm.o.nic  of  electro.n-cyclotro.n 
and  electrostatic  ion-cyclotro-n  oscillations  by  m.eans  of  a m.onochromatic 


f ■>  I • , 


elect rom.aii.-ietic  wave  at  jb  =9.  23  GKz  has  been  observed  exoe rim.entallv. 

''  / o * ' S... 

The  experim.e.nt  has  bee.n  conducted  in  a mdcrowave  sustai.ned  beam,  plasm.a. 
By  switching  the  purnip  power  betwee.n  two  levels,  the  time  evolutic.n  of  the 
excited  electrostatic  ion-cyclotron  wave  is  studied.  In  order  to  explai.n  the 
observed  nonli.ncar  saturated  state  of  the  i.nstability,  the  .ncnli.near  dam.pi.ng 
rate  due  to  two  possible  mecha.nisms  has  been  examined.  Harm.o.nic  gener- 
ation gives  very  small  contributio.ns . therefore,  a.nomalous  diffusion  caused 
by  the  coherent  oscillation  of  the  electrostatic  ion-cyclotron  wave  plays  an 

imiportant  role.  Working  in  the  guiding  center  frame,  the  effective  non- 
linear damping  frequency  pertaining  to  the  diffusion  process  is  obtai.ned. 


